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We obtain constraints on the slope of a universal stellar initial mass function (IMF) over a range 
of cosmic star-formation histories (SFH) using z « OT luminosity densities in the range from 0.2 /im 
to 2.2 /im. The age-IMF degeneracy of integrated spectra of stellar populations can be broken for 
the Universe as a whole by using direct measurements of (relative) cosmic SFH from high-redshift 
observations. These have only marginal dependence on uncertainties in the IMF, whereas, fitting to local 
luminosity densities depends strongly on both cosmic SFH and the IMF. We fit to these measurements 
using population synthesis and find the best-fit IMF power-law slope to be F = 1.15 ± 0.2 (assuming 
dTV/dlogm oc m.-^ for O.5-12OM0 and m~"-^ for 0.1-0.5Mq). This M > 0.5Mq slope is in good 
agreement with the Salpeter IMF slope (F = 1.35). A strong upper limit of F < 1.7 is obtained which 
effectively rules out the Scalo IMF due to its too low fraction of high-mass stars. This upper limit is at 
the 99.7% confidence level if we assume a closed-box chemical evolution scenario and 95% if we assume 
constant solar metallicity. Fitting to the Ha line luminosity density, we obtain a best-fit IMF slope in 
good agreement with that derived from broadband measurements. 

Marginalizing over cosmic SFH and IMF slope, we obtain (95% conf. ranges): figtars — 1.1-2.0 xlO^"^ 
for the stellar mass density; psFR — 0.7-4.1 xlO~^ /i Mq yr^^ Mpc~^ for the star-formation rate 
density, and; — 1.2-1.7 xlO'^^ /iWMpc"^ for the bolometric, attenuated, stellar, luminosity density 
(0.09-5 /im). Comparing this total stellar emission with an estimate of the total dust emission implies a 
relatively modest average attenuation in the UV (< 1 magnitude at 0.2/im). 

Subject headings: stars: formation — stars: luminosity function, mass function — dust, extinction — 
galaxies: luminosity function, mass function — galaxies: stellar content — cosmology: 
miscellaneous 



1. INTRODUCTION 

The luminosity density of the Universe in the range 
from about 0.2 /im to 2.2 /im, mid-UV to near-IR, is dom- 
inated by stellar emission. This wavelength range spans 
the peak in the spectra (/^) of stars with effective tem- 
peratures from about 20000 K down to 2000 K. Therefore, 
measurements of the luminosity density in various broad- 
bands across this range provide a powerful constraint on 
cosmic star- formation history (SFH) and/or a universal 
stellar initial mass function (IMF). 

The stellar IMF describes the relative probability of 
stars of different masses forming (see Gilmore & Howell 
1998, for recent reviews and analyses). Its importance 
crosses many fields of astronomy from, for example, star 
formation (testing theoretical models) to cosmic chemical 
evolution (heavy metal production from high mass stars). 
It is widely used in the study of the SFH of galaxies from 
their integrated spectra. Generally, an IMF is assumed 
and used as an input to evolutionary, stellar population, 
synthesis models, and these models are fitted to integrated 
spectra. 

The first calculation of an IMF was made by Salpeter 
(1955) based on the observed luminosity function of 
solar-neighborhood stars, converting to mass, correcting 
for main-sequence lifetimes and assuming that the star- 
formation rate (SFR) has been constant for the last 5 Gyr. 
Despite the uncertainties in mass-to-light ratios, stellar 
lifetimes and the SFR, this result (a power-law slope of 

^ The Salpeter IMF power-law slope is —1.35 wrt. logarithmic mass bins and —2.35 wrt. linear mass bins. 
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— 1.35 measured from about 0.3 to I5M0) is still com- 
monly used today. ^ 

Measurements of the solar neighborhood IMF were re- 
viewed by Scalo (1986) producing an IMF with a mass 
fraction peak around 0.5-1 Mq (Fig. 1). When applied to 
galaxy populations, this IMF is unable to reproduce Ha 
luminosities (Kennicutt, Tamblyn, & Congdon 1994), and 
applied to cosmic evolution, is unable to match observed 
mean galaxy colors (Madau, Pozzetti, & Dickinson 1998) 
due to a too low fraction of high-mass stars (M > 10 M0). 

In a more recent review by Scalo (1998), he con- 
cluded that the field star IMF was of questionable use 
for a number of reasons (e.g. the derived IMF in the 
range O.9-1.4M0 depends strongly on the assumed solar- 
neighborhood SFH). Instead, he summarized the results 
from studying star clusters in a triple-index power-law 
IMF as an estimate of an average IMF (Fig. 1). He also 
noted that "if the existing empirical estimates of the IMF 
are taken at face value, they present strong evidence for 
variations, and these variations do not seem to depend 
systematically on physical variables such as metallicity 
or stellar density" . The observed IMF variations from 
stellar counts could be largely due to statistical fluctua- 
tions and/or observational biases such as mass segregation 
within a star cluster (Elmegreen 1999). These points, if 
correct, mean that the concept of a universal IMF is highly 
useful in the studies of integrated spectra of galaxies (since 
they are mostly the result of many star-formation regions 
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Fig. 1. — Stellar initial mass functions: mass fraction (per logarithmic mass bin) versus mass for the Salpeter (1955) IMF, the Miller & Scalo 
(1979) IMF, the Scalo (1986, 1998) IMFs, the Kroupa (2001, cqns. 2 and 6) IMFs (left panel) and the IMF parameterization of Equation 3 
(right panel). All are assumed to be valid over the range 0.1-120 Mq. The integral of each curve is set to unity. Note that the slopes of the 
lines are equivalent to 1 — F [= dlog(mniog^) / dlog(m), see Eqn. 3 for definitions]. The F values for M > IM© for the published IMFs are 
1.35 (S55), 1.5/2.3 (MS79), ~2.05/1.5 (S86), 1.7/1.3 (S98) and 1.3 (KOI) (see also Table 2). 



and episodes) but the IMF is still uncertain at the level of 
0.5 in the power-law slopes. By 'universal IMF', we mean 
an IMF that represents the average in a significant ma- 
jority of galaxies and over a significant majority of cosmic 
time. 

Measurements of the cosmic luminosity densities repre- 
sent the A/ A A ^ 6 components of the "Cosmic Spectrum" 
which is the luminosity-weighted, average galaxy spectrum 
(Baldry et al. 2002). Thus, if there is a universal IMF, it 
should certainly apply to this spectrum and constraints 
on the IMF can be obtained by fitting population syn- 
thesis models. These constraints are degenerate with the 
assumed cosmic SFH. However, cosmic SFH is generally 
more accurately quantified than the SFH of any individ- 
ual galaxy other than in the Local Group. For example, 
star-forming galaxies have been observed as far back as 
2; ~ 5 (Dey et al. 1998), and there is evidence for reion- 
ization as early as 2; ~ 20 (Kogut et al. 2003), which for a 
Universe age of 13.7 ± 0.2 Gyr (Spergel et al. 2003) gives a 
time of 12.3-13.7 Gyr since the onset of a significant rate of 
star formation in the Universe. The results for individual 
galaxies arc often limited by age-metallicity degeneracy 
(Worthey 1994) and recent bursts of star formation which 
can disguise their underlying age (Barbaro & Poggianti 
1997). Given our knowledge of cosmic SFH, we can then 
place constraints on a universal IMF by limiting the SFH 
we consider. 

The primary knowledge of cosmic SFH comes from mea- 
suring the comoving density of SFR indicators at various 
redshifts (Madau et al. 1996). These SFR indicators in- 
clude UV luminosities and emission-line luminosities that 
are dominated by light from short-lived high-mass stars. 
The conversion to SFR depends on the IMF. However, the 
relative cosmic SFH is well defined if the same indicator is 
used at each redshift regardless of the assumed IMF. Even 
if the indicator varies (e.g. 0.2 to 0.3 /xm rest-frame UV), 



the derived cosmic SFH is less sensitive to the slope of the 
IMF than the local luminosity densities (0.2 to 2.2/xm). 
It is this principal that enables constraints on a universal 
IMF. 

An analysis of this type was applied by Madau et al. 
(1998) using various luminosity densities (0.15 to 2.2/im) 
spread over a range of redshifts (0-4). They fitted cosmic 
SFHs for three different IMFs. Here, we use a more quan- 
titative approach to constraining a universal IMF slope 
and use more accurate, recent, local {z « 0.1) luminosity 
density measurements. Note that wc assume that there 
is a universal IMF and do not constrain any variation in 
the IMF between different galaxies (see e.g. Wyse 1997; 
Kroupa 2002, for evidence for an invariant IMF). Even if 
there is some variation, the results presented here could be 
regarded as constraints on a luminosity-weighted, average 
IMF. 

To summarize, if the cosmic SFH is assumed to be 

known, based on measurements of various SFR indicators 
with redshift, then local luminosity density measurements 
provide a constraint on a universal stellar IMF. The other 
principal factors to consider are chemical evolution (metal- 
licity) and dust attenuation. 

The plan of the paper is as follows. In Section 2, we de- 
scribe details of recent luminosity density measurements. 
The measurements arc illustrated in Figure 2 and summa- 
rized in Table 1. In Section 3, we describe our modeling 
and fitting procedure. A summary of the parameters used 
in the modeling is given in Section 3.5. In Sections 4 and 5, 
we present our results and conclusions. 

2. LUMINOSITY DENSITY MEASUREMENTS 

In this section, wc summarize the details of various lu- 
minosity density measurements and their conversion to 
AB magnitudes (Oke & Gunn 1983) per comoving Mpc^. 
The zero point of the AB absolute magnitude scale is 



Constraints on an IMF from Luminosity Densities 



3 



-IE 



o 

Oh 



-17 



-16 



-15 



-14 



-13 



TTT 



O Sullivan et a 

□ Madgwick el 
A Norberg et a 
X Blanton et a 
^ Blanton et a" 
A Cole et al. '0 
X Kochanck ct 

□ HiuuiH et al. 



1. '00 - 
al. '02 
1. '02 
I. '01 
I. '03 - 
1 

al. '01 

03 



0.15 0.2 0.3 0.5 0.7 1.0 1.5 
rest -wavelength (/im) 



2.0 



Fig. 2. — Local luminosity densities: absolute AB magnitudes per Mpc3 [(/i.Q^o.Qao) = (1.0,0.3,0.7)] versus wavelength. The vertical 
bars represent uncertainties while the horizontal bars represent FWHM of the bandpasses. See Table 1 and Section 2 for details. We use the 
Sullivan et al. (2000), the Blanton et al. (2003b) and the J- X- band results for our fitting. 



Table 1 

Local luminosity densities from various surveys {0 < z < 0.2) 



band 


reference 


Aeff 


j + 2.5 log ft. 


notes 






(/zm) 


"(AB magnitudes) 




"FOCA 0.2 /xm 


Sullivan ct al. (2000) 


0.20 


-13.83 ±0.25 


z « 0.15, uncorrected for dust 


'^APM 6j 


Madgwick et al. (2002) 


0.46 


-15.53 ±0.10 


z « 0.10 




Norberg et al. (2002) 


0.46 


-15.43 ±0.10 


evolution corrections to ^; = 0.0 


''SDSS u 


Blanton et al. (2001) 


0.35 


-15.25 ±0.25 


z « 0.05 


SDSS g 


55 


0.47 


-16.01 ±0.15 


z « 0.10 


SDSS r 


55 


0.62 


-16.40 ±0.10 


)) 


SDSS i 


55 


0.75 


-16.73 ±0.10 


V 


SDSS z 


55 


0.89 


-17.02 ±0.10 


11 


SDSS o 


Blanton et al. (2003b) 


0.32 


-14.10 ±0.15 


evolution corrections to z = 0.1 


SDSS °-^g 


55 


0.42 


-15.18 ±0.05 


11 


SDSS °-V 


55 


0.56 


-15.90 ±0.05 


11 


SDSS 


55 


0.68 


-16.24 ±0.05 


11 


SDSS 0-1^ 


55 


0.81 


-16.56 ±0.05 


11 


^2MASS J 


Cole et al. (2001) 


1.24 


-16.58 ±0.15 


evolution corrections to z = 0.0 


2MASS Ks 


55 


2.16 


-16.82 ±0.15 


evolution corrections to z = 0.0 




Kochanek et al. (2001) 


2.16 


-16.94 ±0.15 


z ^ 0.03 


Hawaii Kg 


Huang et al. (2003) 


2.16 


-17.57±0.15 


z^0.15 



"luminosity density conversion: j = — 2.51og(Lj// WHz~^ Mpc^"^) + 34.10 
''FOCA magnitudes: assumed toab = "?,2ooo + 2.25 (Milliard et al. 1992) 
''APM magnitudes: assumed toab = bj — 0.08 (see Section 2.2) 

''SDSS magnitudes: assumed toab = {u,g,r,i,z) + (-0.04,0.035,0.015,0.015,0.00) (Blanton et al. 2003b) 
'^2MASS, Hawaii magnitudes: assumed mAB = J + 0.9, mAB = ± 1.8 (Cutri et al. 2001) 
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4.345 X (3631Jy for apparent magnitudes). In 

general, an estimate of the total luminosity density from 
the galaxy population is obtained by an analytical inte- 
gration of the Schechter (1976) function parameters, given 

by 

j = M,- 2.5 log Tf (a + 2)] + C (1) 

where Ff is the gamma function and C represents correc- 
tions from the magnitude system defining the luminosity 
function to total AB magnitudes. The final luminosity 
densities are quoted for a cosmology where (ft, Q,mo , ) 
= (1.0,0.3,0.7) and h = Fo/lOOkms"! Mpc'i. The sur- 
veys select redshifts to a limiting magnitude in the same 
wavelength as that of the luminosity density measurement. 

2.1. Ultraviolet: 0.2i^m 

An analysis of a mid-UV selected redshift survey is pre- 
sented by Sullivan et al. (2000). The imaging for this 
survey covers about 6deg^ to a depth of about 21 AB 
magnitudes using the FOCA balloon-borne telescope (Mil- 
liard et al. 1992). The filter response approximates a 
Gaussian centered at 2015 A with a full-width half max- 
imum (FWHM) of 188 A. The survey area covers four 
fields chosen with very low milky-way (MW) extinction 
[E{B — V)< 0.015]. Follow up spectroscopy was obtained 
using multi-object, optical spectrographs on the WIYN 
and WHT telescopes (Treyer et al. 1998). 

Sullivan et al. measured the local luminosity density 
using a sample of about 200 galaxies (0 < 2: < 0.4) 
to estimate the luminosity function. The Schechter 
parameters without dust-correction were (M* , 0* , a) = 
(—20.59,0.00955,-1.51). Integrating this function gives 
—16.18, and converting to AB magnitudes gives —13.93 
using mAB = ^2000 + 2.25 (Milliard et al. 1992). We also 
correct the measurement to our default world model of 
(fiTOoj^Ao) = (0.3,0.7) from (1.0,0.0) assuming a redshift 
of 0.15. This gives a correction of about -1-0.1 magnitudes. 
The formal uncertainty from the fitting is 0.13 magnitudes. 
However, there are additional uncertainties: absolute cal- 
ibration (< 0.25 mags), MW extinction (< 0.15 mags), 
conversion to total magnitudes (< 0.1 mags), large-scale 
structure, incompleteness corrections, etc. We will assume 
an additional la uncertainty of 0.2 to be added in quadra- 
ture, so that j = -13.83 ± 0.25. 

The limit of this survey is m2ooo = 18.5 which corre- 
sponds to M* galaxies at z « 0.2. With a steep faint-end 
slope of —1.5, the luminosity- weighted mean redshift is 
around 0.15 corresponding to the redshift position plotted 
by Sullivan et al. This is higher than our fiducial redshift of 
0.10 (see below), giving a higher luminosity for any declin- 
ing cosmic SFR at 2 < 0.5. However, this may be counter- 
acted by the lack of MW-extinction and total-magnitude 
corrections. For simplicity and since we do not want to 
assume a cosmic SFH and IMF, we will use the luminosity 
density as it was measured. Note that a couple of galax- 
ies with obvious active galactic nuclei (AGN) characteris- 
tics were removed from their sample. We will assume the 
measured luminosity density does not have a strong, non- 
stellar, AGN component. Neither Sullivan et al. (2000) or 
Contini et al. (2002) found strong evidence for significant 
AGN contamination based on emission-line flux ratios. 

2.2. Optical: 0.3 jim to 0.9 jim 



The Sloan Digital Sky Survey (SDSS; York et al. 2000; 
Stoughton et al. 2002) has imaged over 2000 deg^ in five 
bandpasses {ugriz) with effective wavelengths from 0.35 
to 0.9 /zm to a depth of about 21-22 AB magnitudes. Fol- 
lowup spectroscopy is also included as part of the SDSS 
for various targeting schemes, of which, the main galaxy 
sample (MGS; Strauss et al. 2002) is appropriate for de- 
termining cosmic luminosity densities. The MGS is a 
magnitude-limited galaxy sample (r < 17.77) with a me- 
dian redshift of 0.10. To form effectively complete sam- 
ples in the other four bands, galaxies were selected to 
{u,g,i,z) < (18.4,17.65,16.9,16.5). 

The first luminosity densities from the MGS were pub- 
hshed by Blanton et al. (2001). However, Wright (2001) 
found that these results overpredicted the near-IR lumi- 
nosity density by a factor of 2.3 (compared to Cole et al. 
2001; Kochanek et al. 2001). Since then, better analysis 
techniques, calibration and more data have allowed the op- 
tical luminosity density measurements to be significantly 
improved. We use the results of Blanton et al. (2003b) 
as the basis for our fitting here. One of their approaches 
was to fc-correct and to evolve-correct to a fiducial redshift 
of 0.10. This reduces systematic uncertainties associated 
with these types of corrections because the median redshift 
of the SDSS is at this mark and the luminosity-weighted 
mean redshift is close to it. Thus, their results are most 
accurate for the shifted bandpasses, designated ^'^u, °'^g, 
o.ij,^ 0.1^ o.i^, (rest-frame bandpasses for galaxies at 
z = 0.1). We use the °'^r band as the fiducial band and 
measure all colors with respect to it when comparing syn- 
thetic magnitudes with luminosity densities. We set a Icr 
uncertainty of 0.05 for the g,r, 2,z band measurements. 
This allows for some mis-calibration since the formal un- 
certainties of Blanton et al. are 0.03/0.02 for these bands. 
Using such small errors is appropriate since it is the rela- 
tive o ig^o i i^o.i ^ o.ij, nieasurements that constrain the 
normalized SFH or IMF. In other words, the errors only 
need to represent the uncertainties in the colors and the 
absolute measurements of the luminosity densities need 
not be accurate to this level. An estimated conversion 
of SDSS to AB magnitudes is given and used by Blanton 
et al. (see also Table 1), and the magnitudes used by them 
are assumed to be close enough to total that a correction 
is not applied. 

The new results of Blanton et al. (2003b) are in good 
agreement with &j luminosity densities determined by 
Madgwick et al. (2002) and Norberg et al. (2002). These 
analyses were based on the Automated Plate Measuring 
(APM; Maddox et al. 1990) galaxy catalog with redshifts 
from the 2dF Galaxy Redshift Survey (2dFGRS; Colless 
et al. 2001). The conversion to AB magnitudes of —0.08 
was based on an integration of the 6j curve (P. C. Hewett & 
S. J. Warren 1998 private communication) through a spec- 
trum of Vega (Lejeune, Cuisinier, & Buser 1997, computed 
by R. L. Kurucz). The 6j magnitudes were calibrated to 
total by comparison with deeper, CCD photometry. 

2.3. Infrared: l.Ofim to 2.5 jim 

The Two Micron All Sky Survey (2MASS; Skrutskie 
et al. 1997) has imaged the whole sky in the J, H and 
Kg bands to a depth of 15-16 AB magnitudes (14.7, 13.9, 
13.1 Vega mags, respectively). Cole et al. (2001) matched 
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the second incremental data release to redshifts obtained 
by the 2dFGRS. With this data set, they determined the 
local J and K band luminosity functions. With k and 
evolution corrections to ^; = 0, the Schechter parameters 
were (-22.36,0.0104,-0.93) and (-23.44,0.0108,-0.96) 
for the J and K bands, respectively. Integrating these 
functions gives —17.36 and —18.50. The conversions to AB 
magnitudes are taken as ruAB = J+0.9 and K + 1.8 (Cutri 
et al. 2001) and the conversion to total magnitudes is es- 
timated to be between —0.08 and —0.15 (we take —0.12). 
The uncertainties come from Poisson noise, the absolute 
calibration, the conversion to total magnitudes and large- 
scale structure. These are not all well defined but we will 
be conservative and use 0.15 for the la uncertainty so that 
j = -16.58 ±0.15 and j = -16.82 ±0.15 for the J and K 
bands, respectively. 

The /C-band luminosity function was also determined 
by Kochanek et al. (2001) using 2MASS imaging. Here, 
they determined the luminosity density using a shal- 
lower sample but with greater sky coverage for the red- 
shifts. Their best estimate was obtained by summing sep- 
arate luminosity functions for late and early-type galaxies: 
(-22.98,0.0101,-0.87) and (-23.53,0.0045,-0.92). The 
total luminosity density is then —16.94 in AB magnitudes 
after appropriate corrections (to AB, as above; to total, 
—0.21). This is in good agreement with the Cole et al. 
result. 

A deeper /iT-band survey covering about 8 deg^ was re- 
cently analysed by Huang et al. (2003). Imaging was taken 
with the University of Hawaii telescopes at Mauna Kea 
Observatory (Huang et al. 1997) with redshifts obtained 
using the 2dF facility on the A AT. The best- fit Schechter 
parameters were (—23.70, 0.0130, —1.37), giving a luminos- 
ity density of j = —17.57 after correcting to AB magni- 
tudes. Even with a conservative error of 0.15, this result 
is > 3(7 discrepant from the result of Cole et al. (2001). 
This discrepancy of 0.75 magnitudes is too large to be 
explained by cosmic evolution since the if -band luminos- 
ity density is dominated by the older stellar populations. 
Analysis by Huang et al. (2003) suggests it could be due 
to a local underdensity. We note that ncar-IR luminosity 
densities have larger uncertainties due to large-scale struc- 
ture than bluer measurements. In underdense regions, 
the SFR per galaxy is higher (as noted by color-density 
and morphology-density relationships, e.g. Blanton et al. 
2003a) which counteracts the effect on the UV luminosity 
density. 

Given the above discrepancy in the near-IR luminosity 
densities, firstly, we will consider the Cole et al. results 
separately from the Huang et al. results in our fitting, and 
secondly, we will use an average iv'-band luminosity. 

3. MODELING 

We fit the data with synthetic magnitudes calculated 
using the PEGASE.2 evolutionary synthesis code (Fioc & 
Rocca-Volmerange 1997, 1999) and integrating the spectra 
through the filter response curves (Fig. 3). The responses 
for the SDSS are taken from Stoughton et al. (2002), the 
FOCA filter curve is assumed to be Gaussian (Milliard 
et al. 1992) and the near-IR filter curves are taken from 
Cutri et al. (2001). We run the PEGASE models with 
nebular emission but 'no extinction' except we modified 



the code so that the absorption of Lyman continuum pho- 
tons by dust was still included (following the prescriptions 
of Spitzer 1978). This was done because our dust atten- 
uation parameterization (Sec. 3.4) does not apply to this 
'pre-extinction' of nebular continuum and line emission. 
The prescription for the models is described below includ- 
ing SFH, IMF and chemical evolution. We include further 
effects of dust attenuation on the output spectra. 

3.1. Cosmic SFH 

For our modeling of cosmic star-formation history we 
use a 'double power-law' parameterization (Baldry et al. 
2002): 

(l + z)^ forO<z<l 
(1 + z)" for 1 < .2 < 5 



SFR oc 



(2) 



The SFRs from the two power-laws are the same at ^; = 1 
and star formation is started at z = 5. Clearly, there are 
many other possible cosmic SFHs other than defined by 
these two parameters. However, the resulting synthetic 
spectra are highly degenerate with SFH. We chose this 
two-parameter model because the /? parameter provides a 
good match to direct measures of SFR at z < 1 and the a 
parameter allows us to add high-redshift star formation in 
a well-defined way. We assume a cosmology correspond- 
ing to {h, il„j„, ilAo)= (0.7,0.3,0.7). Most importantly this 
determines the timescale for cosmic SFH measurements. 

Examples of this parameterized cosmic SFH are shown 
in Figure 4 with a timescale for the h = 0.7 cosmol- 
ogy. Note that direct measures of the relative cosmic 
SFH with redshift do not depend on Hq. The Hubble 
constant changes the scaling factor at all redshifts by the 
same amount. However, our modeling of the cosmic SFH, 
to fit to local luminosity densities (fossil cosmology), does 
depend on Hq because of the timescale dependence. The 
synthetic spectra are calculated at 11 Gyr corresponding to 
the fiducial redshift, z = 0.1, with the chosen cosmology 
and Zform = 5. 

There is overwhelming evidence for a rise in SFR to 
z = 1 from a variety of indicators (e.g. Haarsma et al. 2000; 
Hammer et al. 1997; Lilly et al. 1996; Rowan-Robinson 
et al. 1997). Recent estimates from compilations of mea- 
sures of luminosity density with redshift out to 2 = 1 give, 
for example, /3 = 2.7 ± 0.7 (Hogg 2002). If we take the 
formal 2(7 range from Hogg, we obtain a range of 1.3~ 
4.1. However, some of the measurements compiled by him 
could be biased toward a steeper evolution (e.g. Lilly et al.) 
because of selection effects. Cowie, Songaila, & Barger 
(1999) found a shallower evolution of /3 = 1.5 ± 0.5 based 
on rest-frame UV selection at all redshifts. To encompass 
the "evidence for a gradual decline" (Cowie et al.) to a 
steep decline (Lilly et al.), we will consider /3 to be in the 
range 0.5-4.0. 

High redshift (z > 1) star formation is parameterized by 
a and zibnn = 5. We do not try to exactly match high red- 
shift measurements. The uncertainties are still large due 
to, for example, dust and surface brightness corrections. 
By taking a from —2 to 2, we can approximate the effect 
of high-redshift star formation from a rapid decline at high 
redshift (Madau et al. 1996) to a flattening (Steidel et al. 
1999) to a significant rise (Lanzetta et al. 2002). Note that 
the effect of any significant star formation in the ~1 Gyr 
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Fig. 3.— FOCA, SDSS and 2MASS filters. The SDSS filters are shown 'shifted' to z = 0.1 since we are fitting to the results of Blanton 
et al. (2003b) . Each curve is normalized so that the integral of the normalized transmission T d In A is equal to unity. Thus the height of each 
curve is related to its resolving power (A/ A A). 
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Fig. 4. — Examples of the parameterized cosmic SFH. The SFR is proportional to (1 + z)'^ for z < 1 and (1 + z)" for 1 < x < 5. Each 
curve is normalized so that the total star formation between z = 5 and z = is unity. The timescale is for (ft., Hmo, S^Ao)= (0.7,0.3,0.7). Six 
SFHs are shown with two values of /3 (1,3) for each of three values of a (—2,0,2). 



between z = 20 and z = 5 can be approximated by an 
increase in a. 

3.2. Universal IMF 

Significant degeneracies exist with modeling the spectra 
of galaxies, e.g. age-metallicity. Given these uncertainties, 
modelers have generally assumed an IMF, typically the 
Salpeter (1955) IMF for modeling galaxy spectra. Here, 
we will assume a universal IMF (constant with time and 
environment) but we will consider different IMF power- 
law slopes. A double power-law IMF is chosen based on 
the rationale that there is a clear change in slope around 
0.5 Mq but no definitive change at higher masses (Kroupa 
2001). Our parameterization is^ 

m~°-^ for 0.1 < m < 0.5 , . 

niogm oc 0.5 < m < 120 ^"^^ 

where niog m d log m is the number of stars with logarithm 
of the mass in the range log m to log m + d log m, m is in 

^ An equivalent formalism for our upper-mass slope is dN/d log m oc 
follows the nomenclature outlined by Kennicutt (1998). 



units of solar masses with 0.1 M© and 120 M© being the 
IMF cutoffs. For the Salpeter (1955) IMF, F = 1.35 in this 
formalism except, traditionally, the slope continues down 
to 0.1 M0. A comparison between different IMFs is shown 
in Figure 1. Note that despite the average IMF of Scalo 
(1998) using a break at 1 and 10 M©, figure 5 of that pa- 
per appears consistent with a single break at 0.5-0.8 M©. 
Thus, as with Kroupa (2001), there is no strong evidence 
to rule out a double-power law IMF with a break at 0.5 Mq 
being an adequate approximation of a universal IMF. For 
the low-mass slope, we use the average of the two Kroupa 
equations. We do not vary this slope since our results are 
less sensitive to this than the upper-mass slope. 

3.3. Chemical evolution 

For chemical evolution, we incorporate the closed-box 
evolutionary model within PEGASE which uses the mod- 
els of Woosley & Weaver (1995) to estimate metallicity 
production via supernovae. The metallicity is controlled 

where N is the cumulative number of stars up to mass m. This 
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by a parameter r (Baldry et al. 2002) that represents the 
total maws of stars formed between z = 5 and z = divided 
by the total amount of gas initially available. Higher val- 
ues of r produce higher metallicity. The parameter can be 
greater than unity because of recycling of material. Since 
we are comparing different IMFs, we do not quote r val- 
ues but rather the value for the metallicity [Z) averaged 
on the luminosity at z = 0.1 (the Z-r dependence varies 
with IMF). Solar metallicity is considered to be Z = 0.02 
in the PEGASE models. 

3.4. Dust attenuation 

We wish to estimate the effective dust attenuation for 

the cosmic spectrum covering the UV to near-IR. How- 
ever, we can not consider a standard slab or screen model 
as we would for an individual galaxy because the contri- 
butions to the luminosity densities come from many types 
of galaxies. Kochanek et al. (2001) estimated 54% and 
46% contributions to the JT-band from late- and early-type 
galaxies, respectively. The ratio is similar in the visible 
red bands from the SDSS survey though it depends on the 
chosen dividing line in color or morphology (Blanton et al. 
2003a) . From Madgwick et al. (2002), 61% and 39% are 
the fractional contributions to the 6j band from late and 
early types, respectively, based on spectroscopic classifica- 
tion (assuming Types 2-4 for late. Type 1 for early). Late- 
type galaxies (Sa-Sd and starbursts) contribute about 90- 
95% of the light in the 0.2 jim UV (estimated from table 7 
of Sullivan et al. 2000). This is consistent with the re- 
sults obtained by Wolf et al. (2003) for the 0.28 /xm UV 
at z = 0.3, with ~ 75/90% contribution from late spectral 
types (3-4/2-4). Using the inclination- averaged attenua- 
tions of various galaxy types given by Calzetti (2001, ta- 
ble 3) and averaging over suitable distributions at each 
wavelength, we obtain cosmic spectrum attenuations in 
magnitudes of 1.1-1.35 at 0.15 /itm, 0.4-0.55 at 0.45 /xm, 
0.2-0.3 at 0.8 ^m and 0.05 0.1 at 2.2 /im. 

To incorporate dust attenuation, we use a power law. 
The attenuation in magnitudes as a function of wavelength 
is approximated by 

where Ay is the attenuation at the fiducial wavelength of 
\n. We use 0.56 /im since it matches the effective wave- 
length of our fiducial '^■^r band and it is close to the stan- 
dard V band (0.55 /zm). 

If we fit to the cosmic spectrum attenuations, estimated 
above, we obtain {Ay,n) « (0.35, 1.0) for our fiducial dust 
model. This curve is within the ranges in attenuation es- 
timated above at each wavelength. Chariot & Fall (2000) 
found that for an average star-burst galaxy, n was ap- 
proximately 0.7 in a UV-to-visible attenuation law. Their 
model also included increased effective absorption in birth 
clouds which had a finite lifetime. Here, we assume that 
the the average of all galaxies (the luminosity densities 
considered here) can be approximated by a single effective 
optical depth. The average spectrum of the Universe is 
not that of a star-burst galaxy and the broadband colors 
are not strongly affected by emission lines. The luminosity 
density attenuations are naturally steeper than n = 0.7 be- 
cause of the increasing contribution of less dusty ellipticals 



as the wavelength is increased. Rather than restricting our 
fitting to a single value of n, we allow for a range from 0.8 
to 1.2. 

For Ay , we allow for a range from 0.2 to 0.55 magnitudes. 
For n = 1.0, this is equivalent to a range in A2000 from 0.55 
to 1.55 magnitudes. This encompasses the difference be- 
tween the uncorrected and dust-corrected luminosity den- 
sities of Sullivan et al. which amounts to 1.3 magnitudes. 
To avoid excessive attenuation at UV and near-IR wave- 
lengths, we also set 0.5 < A2000 <i 1-6 and Ak < 0.15 
which reduces the A^ range away from n = 1.0. In gen- 
eral, we marginalize over attenuation, i.e., we chose the 
best fit parameters within the defined ranges that mini- 
mize when fitting the SFH and IMF. 

3.5. Summary of parameters 

In this section, we summarize the important parameters 
and ranges considered in our analysis and show the effect 
of varying some of the parameters. 

• Synthetic spectra. PEGASE.2 (Fioc & Rocca- 
Volmerange 1997, 1999): Padova tracks (Bressan 
et al. 1993), spectral libraries of Lejeune et al. 
(1997) and Clegg & Middlemass (1987). The Leje- 
une et al. library is principally derived from Kurucz 
(1992) model atmospheres. Nebular continuum and 
line emission is also included. 

• Cosmology. {h.^moi^Ko) = (0.7,0.3,0.7) (except 
luminosity densities are quoted for /i = 1, or, equiv- 
alently as j + 2.5 log /i). These values are fixed in 
our analysis. 

• Cosmic SFH (Sec. 3.1, Eqn. 2). [3 between [0.5,4.0], 
a between [—2,2], Zform = 5. Low-redshift {z < 1) 
star formation is parameterized by SFR a {1 + z)^ 
and high-redshift by (1 -|- z)". 

• Universal IMF (Sec. 3.2, Eqn. 3). Power-law 

slopes of —0.5 (for m < 0.5 M0) and — F between 

[-0.8, -1.8] (for TO > 0.5 Mq) with (mmin, Wmax) = 

(0.1,120). The power-laws are with respect to logm 
so that the Salpeter IMF slope is —1.35. 

• Chemical evolution (Sec. 3.3). Closed-box approxi- 
mation with Z between [0.008,0.05] (Z0==O.O2) from 
r between [0.1,1.4]. The final luminosity- weighted 
metallicity Z is principally a function of r and the 
IMF. 

• Dust attenuation (Sec. 3.4, Eqn. 4). Power-law 
slope of —n between [—0.8, —1.2] with Ay between 
[0.2,0.55]. The ranges are also constrained by 0.5 < 
^2000 < 1-6 and Ak < 0.15. 

Figure 5 shows the effects of varying one parameter on 
the synthetic spectra (normalized at the °'^r band) with re- 
spect to a model of (F, f3, a, Z, Ay, n) = (1.3, 2.5, 0, 0.018, 
0.35, 1.0). Over our chosen parameter ranges, varying the 
IMF slope r has the largest effect on the UV — "'^r colors, 
with (} having the second largest effect. The metallicity 
has the largest effect on the near-IR — °-^r colors. 
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Fig. 5. — Tile effect of varying selected parameters individually: absolute magnitude versus wavelength. The lines represent synthetic 
spectra smoothed to A/AA ~ 10. The dotted line is derived from the same set of parameters in every panel [(F, /3, a, Z, A^, n) = (1.3, 2.5, 0, 
0.018, 0.35, 1.0)]. Ail the spectra are normalized to —15.90 in the "'^r band filter (Blanton et al. 2003b). The symbols represent luminosity 
density measurements; the diamonds for Sullivan et ai. and Blanton et al.; the triangles for Cole et al., and; the square for Huang et al. 



4. RESULTS 

First, we fit to the UV-to-optical luminosity density 
measurements of Sullivan et al. (2000); Blanton et al. 
(2003b) with different near-IR measurements, considered 
separately: (a) using the Cole et al. (2001) i^-band result 
(compilation designated SBCx ); (b) using the Huang et al. 
(2003) result (compilation SBH), and; (c) using the Cole 
et al. J-band result (compilation SBC./) and not includ- 
ing the z-band measurement which is included in the first 
two compilations. Examples of fitted spectra are shown in 
Figure 6. This shows that the near-IR luminosity density 
measurements can be fitted, primarily by varying metal- 
licity, while the UV-to-optical spectrum remains approx- 
imately the same. Within the range of metallicity con- 
sidered here (0.4 Z0-2.5 Z©), the best fit is obtained to 
the SBCa: data. Here, all the measurements are within 
about 1(7 of the synthetic magnitudes. The new results of 
Blanton et al. (2003b) and this analysis resolves the major 



discrepancy noted by Wright (2001) between the optical 
and near-IR luminosity densities. However, there is a sig- 
nificant discrepancy between the Hawaii and the 2MASS 
X-band luminosity densities and a discrepancy between 
the J-band result and the z-band/if-band results. 

Neither of the infrared surveys is ideal for comparison 
with the SDSS survey results. The 2MASS survey is too 
shallow (median z ~ 0.05 at the limit of the extended 
source catalog) for direct comparison with the 2; = 0.1 
analysis and the Hawaii survey goes deeper but has a sig- 
nificantly smaller area (8deg2). Note that the Sullivan 
et al. result only uses about 200 galaxies but it is con- 
sistent with the °-^w band result. Improved surveys are 
needed to resolve the near-IR luminosity density discrep- 
ancies and to reduce the uncertainties in UV luminosity 
densities. 

We first look at constraints on the metallicity and IMF 
for a given cosmic SFH (marginalizing over dust attenua- 
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Fig. 6. — Example fits to data from SBCk (UV, ugriz and K triangle; solid line), SBH (UV, ugriz and K square; dashed line) and SBCj 
(UV, ugri and J; dotted line). The lines represent the fitted synthetic spectra smoothed to A/AA ~ 10. The parameters for these fits were 
lr,l3,a,Z,Ay,n) = (1.1, 4.0, 0, 0.022, 0.20, 0.9) for SBCic, (1.1, 3.5, 0, 0.048, 0.44, 0.8) for SBH and (1.1, 4.0, 0, 0.009, 0.22, 0.8) for SBCj. 
See text for details. 



tion). We compute sjid determine the confidence levels 
of Ax^ = (1.0, 2.3, 6.2, 11.8) corresponding to (68.3%) 1- 
parameter and (68.3%, 95.4%, 99.73%) 2-parameter confi- 
dence limits. The resulting contours are shown in Figure 7 
for a cosmic SFH with 6 = 2.5 and q = (Hogg 2002; 
Steidel et al. 1999). As expected, the best fit metallicity 
is dependent on the near-IR data: the SBCk data has a 
best fit metallicity around 1-1.5 Z©; SBH data, > 2Zq, 
and; SBCj data, < IZ©. The SBCi^ set of data is in 
agreement with solar neighborhood measurements of the 
metallicity which give an average close to solar (Haywood 
2001). The SBH data is in disagreement with solar metal- 
licity at the 99.7% confidence level. 

For the rest of the paper, we use an average estimate of 
the ii'-band luminosity corresponding to 

j ^ -16.88 ±0.25 (5) 
from Cole et al. (2001) and Kochanek et al. (2001). The 
uncertainty was increased, and the J-band was not in- 
cluded, so that the constraints on the IMF were not de- 
pendent on the discrepancies noted above. This Itr range 
in the K-hand luminosity density is also similar to the 
range determined by Bell et al. (2003). The metallicity- 
IMF contours for this data set (compilation designated 
SBav) are also shown in Figure 7. The best fit IMF slope, 
r = 1.3 ± 0.1, shows the tightness of the constraint if the 
SFH is known. The best-fit metallicity is greater than so- 
lar but for other cosmic SFHs (e.g. /3 = 2.5 and a = 2), 
the best fit is around solar. In general, there is minimal 
degeneracy between F and metallicity for the data, i.e. 
the choice of metallicity does not significantly affect the 
constraints on the IMF slope. 

Note that there is an upper limit on the metallicity as 
a function of F due to the closed-box model (Fig. 7). In- 
sufficient metals can be produced in the available time to 
raise the average metallicity above a certain limit and this 
limit decreases as the fraction of high-mass stars in the 
IMF decreases. This chemical evolution limit is derived 
from the PEGASE code using the upper yields of Woosley 



& Weaver (1995). 

The pattern shown in Figure 7 is similar as /3 is varied 
except the best-fit F shifts. To illustrate this, we now fix 
Z equal to 0.02 for the SBav fit and show joint confidence 
levels in f3 versus F in Figure 8. We plot the result for 
three values of a. All the plots show a broad degeneracy 
with a slope of —7 ± 1. In other words, the best-fit F de- 
creases by about —0.14 for each -|-1 increase in /3. Note 
that the metallicity was chosen to be consistent with the 
solar neighborhood and so that the closed-box model was 
valid out to F of 1.8 with only minimal extrapolation (see 
Fig. 7). 

The choice of high redshift star formation (defined by a) 
makes little difference on the best-fit F (Fig. 8). Marginal- 
izing over SFH {(3 & a), we obtain a best fit for F in the 
range 0.85-1.3 (68% conf.) and a strong upper limit of 
F < 1.7 (99.7% conf.). This is our principal result. In 
other words, assuming a declining SFR from z = 1 to the 
present day (/3 > 0.5) then the present luminosity densi- 
ties, in particular, the UV to optical colors mean that the 
upper- mass IMF slope cannot be steeper than 1.7. The 
caveats are: (i) PEGASE evolutionary tracks and stellar 
spectra arc sufficiently accruate; (ii) cosmic SFH can be 
approximated with the double power-law and the look- 
back time to z = 1 is close to 7.7 Gyr; (ui) there is a close 
to universal IMF with a near unbroken slope above 0.5 M©; 
(iv) chemical evolution can be approximated by a single 
metallicity for each epoch with the metallicity increasing 
with time (closed-box approximation); (v) the average ef- 
fect of dust can be approximated by a power-law within 
the ranges considered; (vi) the Copernican principle, that 
we are in no special place in the Universe, applies; (vii) 
the luminosity density uncertainties can be approximated 
as Gaussian. 

We now turn to look at varying a couple of these as- 
sumptions, related to metallicity and dust, in Sections 4.1 
and 4.2. In Section 4.3, we calculate the stellar mass, 
SFR and bolometric densities from our models and, in 
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Fig. 7. — Joint confidence levels over metallicity and IMF slope assuming a cosmic SFH (/3 = 2.5, a = 0). Left: The contours represent 
confidence levels for the SBCx data (solid lines), SBH data (dashed lines) and SBCj data (dotted lines). The levels correspond to 68% 
1-para. and 95% 2-para. confidence limits, marginalized over dust attenuation. The thick gray line represents the approximate upper limit 
on the luminosity- weighted Z due to chemical evolution (from the closed-box model of PEGASE). Any values above this line are linearly 
extrapolated. Right: The same except the confidence limits are for the SBav data which uses an average near-IR luminosity estimate (Eqn. 5), 
and the levels represent the 68% l-paxa. and 68%, 95%, 99.7% 2-paj:a. confidence limits. 



Section 4.4, we apply constraints using measurements of 
Ha luminosity density. 

4.1. Metallicity approximations 

Wc have used the 'closed-box' approximation for the 
evolutionary synthesis in the above analysis. This is valid 
if the dominant star formation at each epoch is taking 
place in an average chemical environment at that epoch. 
This scenario could result from complete mixing between 
different environments. Hot, young stars have a higher 
metallicity on average than cool, old stars. Another com- 
monly used assumption in evolutionary synthesis is the 
'constant-metallicity' approximation. This scenario could 
result from no mixing between environments. The metal- 
licity evolves independently and rapidly in each separate 
star-forming region with a characteristic or average metal- 
licity representing all epochs (separate regions) of star for- 
mation. 

Since the Universe has neither complete or no mixing, 
something between these two approximations might be ex- 
pected.^ Examples of differences between these are shown 
in Figure 9. Though the effect on the broadband colors de- 
pends on the SFH and IMF, the general effect is to increase 
the UV and near-IR fluxes with respect to the visible fluxes 
for the constant-.^ approximation. This approximation 
leads to lower metallicities for young/hot stars and higher 
metallicities for old/cool stars compared to the closed-box 
evolution scenario (matching Z to constant Z). 

Some results for the constant-metallicity approximation 
are shown in Figures 10 and 11 for the SBav data. If we 
chose solar metallicity then the results are similar to the 
closed-box approximation except F is in the range 1.05- 
1.45 (68% conf.) and F < 1.7 at 95% confidence. The 
latter confidence level holds for Z < 0.025. 

Table 2 shows a comparison between different pub- 
lished IMFs with multiple power-law slopes. We use the 

^ We do not consider infall or outflow models which are often considered 



constant-metallicity approximation for this comparison to 
avoid the additional complication introduced by metal pro- 
duction in a closed-box model. With this comparison, the 
Kroupa, Tout, & Gilmore (1993); Miller & Scalo (1979); 
Scalo (1986) IMFs are strongly rejected. This is consistent 
with their average slopes, over 1-10 M0 and 10-120 M0, 
being F > 1.7. Note that our modeling cannot strongly 
distinguish between different slope changes below IMq. 

4.2. Dust attenuation approximations 

We have used a power law to describe average dust at- 
tenuation based on estimating the distributions of galaxy 
types at each wavelength (Sec. 3.4, A„ = 0.2-0.55, n = 
0.8-1.2). Chariot & Fall (2000) showed that star-forming 
galaxies were consistent with a shallower or grayer slope 
of n = 0.7. If we take the approximate distribution in the 
attenuation parameter measured by Chariot et al. (2002), 
Ay = 0.8 ± 0.3, and we also take n = 0.7 ± 0.1, both 
with normal distributions (cutoff > 0), then the average 
effective parameters over many galaxies are equivalent to 
{Ay,n) ~ (0.75,0.65). In other words, there is a slight re- 
duction in the effective attemiation and a slight flattening 
of the curve because fluxes are averaged rather than mag- 
nitudes. The attenuation for these average parameters at 
0.2 /xm is 1.45 magnitudes, not far from the average atten- 
uation estimated by Sullivan et al. (2000) of 1.3. Their 
attenuation corrections were based on Balmer line mea- 
surements with conversion to UV attenuation using the 
Calzetti (1997) law. 

In Figure 12, we compare this star- forming galaxy at- 
tenuation law (0.75,0.65) with the average attenuation es- 
timated in Section 3.4 (0.35,1.0) and a steeper MW ex- 
tinction law {Ay = 0.5, Pei 1992). The best fit regions for 
these fixed attenuation laws lie primarily within the best 
fit region based on marginalizing over metallicity. Thus, 
our results are not strongly dependent on the assumptions 

for individual galaxies. 
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Fig. 8. — Joint confidence levels over fl/a and IMF slope for three different values of a (—2,0,2) / one value of /9 (2.5) and for a chemical 
evolution closed-box approximation with Z = 0.02. The contours represent 68% l-para. and 68%, 95%, 99.7% 2-para. confidence levels for 
the SBav data, marginalized over dust attenuation. 



about dust and we have chosen a fairly generous range of 
parameters for our principal fitting. Note that the grayer 
law of (0.75,0.65) produces an unrealistically high atten- 
uation in the ii'-band of 0.3 magnitudes and the MW ex- 
tinction law {n ~ 1.5 visible-to-near-IR) is naturally too 
steep because it is a foreground-screen law. In Section 4.3, 
we analyse the consequences of our dust model for IR lu- 
minosity density due to dust emission. 

4.3. Stellar mass, SFR and bolometric densities 

We can calculate some derived physical properties of the 
Universe (at z « 0.1). The mass density in stars, fistars, is 
in the range 1.1-2.0 xlO~^ h~^, and; the SFR density, 
PSFR, is in the range 0.7-4.1 xl0~^ /iMg yr'^ Mpc~^. 
These ranges represent 95% confidence limits marginal- 
ized over IMF and cosmic SFH but restricted to near- 
solar metallicity models {Z ot Z = 0.015-0.025). How- 
ever, the results depend strongly on the low-mass end of 
the stellar IMF which is not constrained by our analy- 
sis. As a test of the varying the low-mass end of the 
IMF, we also computed the ranges using the best-fitting 
IMFs and cosmic SFHs described in Table 2 (with < 95% 



conf. of rejection). From these, Ogtars is in the range 0.8- 
2.5 xlO~^ h~^, and; psFR is in the range 1.1 4.3 xl0~^ 
/i Mq yr^^ Mpc^''. The lower stellar mass densities are 
derived from IMFs with rm<i < 0.5 (Kennicutt, Salpeter 
mod. B, Scalo 1998) while the higher mass densities are 
derived from the Kroupa (2001) B IMF. The uncertainty 
in the current SFR density is not increased compared to 
the measurement based on the Equation 3 parameteriza- 
tion. The lowest SFR densities (< 0.01) only occur in our 
models with F < 1.1 and none of the published IMFs in 
the table have r„i>i < 1.2 which explains the lack of low 
SFR densities based on those IMFs. 

These results are in good agreement with the results, 
based on the Kennicutt IMF, of Cole et al. and Baldry 
et al., and with the results, based on a modified Salpeter 
IMF, of Bell et al. (2003). They are generally not in agree- 
ment with results based on the Salpeter IMF extending 
down to 0.1 M0 because such an IMF produces a high mass 
density from low mass stars that have minimal impact on 
the luminosity densities. This type of universal IMF is 
ruled out by stellar counts in the MW (e.g. Scalo 1986, 
1998) and by analysis of the dynamics of spiral galaxies 
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Fig. 10. — Joint confidence levels over metallicity and IMF slope assuming a cosmic SFH (/3 = 2.5, a = 0) for the constant-metallicity 
approximation. See Figure 8 for contour meanings. 



(e.g. Bell & de Jong 2001). 

The total, bolometric, attenuated, stellar, luminosity 
density (0.09-5 /im) is determined from the models to be 
in the range 1.2-1.7 xlO^^ ftWMpc"^ (95% confidence). 
We can also estimate the total, bolometric, luminosity den- 
sity 5-1000 /im) due to dust absorbing and re-emitting 
stellar light. This depends critically on our dust model. 
Prom the best-fit models after marginalizing over dust- 
model parameters, the total is in the range 0.3-1.5 xlO^^ 
h W Mpc~^ which corresponds to 20-50% of the unatten- 
uated stellar light being absorbed. If we restrict our dust 
model to {Ay,n) = (0.35,1.0), the ranges are 0.55-0.95 
xlO^^ /iWMpc"^ and 30-40%. 

From Saunders et al. (1990), the local, far-IR, 42- 
122 /xm, luminosity density is in the range 0.17-0.26 x 10"^^ 
/iWMpc"^ (±2c7 range). This is significantly lower than 



the energy predicted by our dust model. However, a cor- 
rection to total dust emission needs to be applied. Using 
the infrared energy dust models of Dale & Helou (2002), 
corrections from the 42-122 band emission to the total 
dust emission range from 1.9 to 2.6. If we assume this rep- 
resents the systematic uncertainty in the luminosity den- 
sity correction, then the total bolometric dust emission is 
in the range 0.3 0.7 xlO''''' hMVMpc^^ (scaUng from the 
Saunders et al. result). Our fiducial dust model has a 
range in total bolometric emission that overlaps with this 
estimate. 

This upper limit of 0.7 xlO^^ /iWMpc"^ for the dust 
emission favors models with lower attenuation and/or 
lower luminosity densities around 0.2 /xm. From our fit- 
ting, after marginalizing over dust-model parameters, we 
obtain j42ooo ^ 1- This is marginally inconsistent with 
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r (IMF slope) r (IMF slope) 

Fig. 11. — Joint confidence levels over /3 and IMF slope for two different values of a (0,2) and for a constant-metallicity approximation with 
Z = 0.02. See Figure 8 for contour meanings. 



Table 2 

Comparison between published IMFs 



reference 


"power-law slopes (— F) 
0.1-0.5 M0 0.5-1 M0 


in four 
1-10 M 


mass ranges 
Q 10-120 Mq 


''confidence of 
rejection 


best-fit SFH (3 
(for a = 0, 2) 


this paper 


-0.5 


-1.2 


-1.2 


-1.2 


< 68% 


3.5, 2.5 


Kennicutt (1983) 


-0.4 


-0.4 


-1.5 


-1.5 


80% 


2.0, 1.0 


Kroupa (2001) A 


-0.3 


-1.3 


-1.3 


-1.3 


<68% 


3.0, 2.0 


Kroupa (2001) B 


-0.8 


-1.7 


-1.3 


-1.3 


<68% 


3.0, 2.0 


Kroupa et al. (1993) 


-0.3 


-1.2 


-1.7 


-1.7 


98% 


1.5, 0.5 


Miller & Scab (1979) 


-0.4 


-0.4 


-1.5 


-2.3 


98% 


1.5, 0.5 


Salpeter modified A 


-0.5 


-1.35 


-1.35 


-1.35 


<68% 


3.0, 1.5 


Salpeter modified B 


-0.5 


-0.5 


-1.35 


-1.35 


<68% 


3.0, 2.0 


'^Scalo (1986) 


-0.15 


-1.1 


-2.05 


-1.5 


99.9% 


0.5, 0.5 


Scale (1998) 


-0.2 


-0.2 


-1.7 


-1.3 


90% 


1.5, 0.5 



"All the IMFs are assumed to be valid from 0.1 Mq to 120 Mq. 

''The IMFs were compared by marginalizing over 24 SFHs (/? = 0.5-4.0, step 0.5, for a = —2,0,2) using the constant- 
metallicity approximation with Z = 0.02. The confidence of rejection is with respect to the best-fit IMF in this table. 
■^The power-law slopes shown for the Scalo (1986) IMF are an approximation from a fit to the mass fractions (Fig. 1). 



the Sullivan et al. (2000) estimate of the effective attenu- 
ation, 1.3, which is why we have not used the estimated 
total dust emission to constrain our dust model. This dis- 
crepancy could be resolved if the UV luminosity density 
and/or attenuation were overestimated,'* and/or the to- 
tal IR plus sub-mm flux was underestimated perhaps due 
to a population of galaxies with colder dust than those 
detected by 60 /xm surveys. Note also that Buat & Bur- 
garella (1998) found an average attenuation of 1.2 but this 
may not be inconsistent with our attenuation limit since it 
represents a limit on the luminosity-weighted average by 
flux. Future analyses could use IR and sub-mm luminosity 
density measurements to better constrain a dust model. 

4.4. Ha luminosity density 



The Ha nebular emission comes from reprocessed Ly- 
man continuum photons. Therefore it provides a measure 
of the UV flux blueward of 0.1 /um. Here, we consider mea- 
surements of the attenuation-corrected Ha luminosity den- 
sity, relative to the '^'^r band, in comparison with model 
predictions. This emission-line attenuation is significantly 
higher than for the stellar light at the same wavelength 
(Calzetti, Kinney, & Storchi-Bergmann 1994) but it can 
be estimated using the Balmer decrement (Ha/H/? « 2.85 
for case B recombination, Hummer & Storey 1987). 

Three attenuation-corrected Ha luminosity density 
measurements are summarized in Table 3. We do not 
quote error bars as the uncertainties are dominated by 
systematics. These include (i) subtraction of the stel- 
lar contamination which, in particular, affects the mea- 
surement of H/3 and thus the estimate of the attenuation 



* However, recent analysis of the FOCA redshift survey with new redshifts and fc-corrections gives a slightly more luminous 0.2 /nm luminosity 
density by about 0.2 magnitudes (M. Sullivan 2003 private communication). 
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r (IMF slope) 



Fk;. 12. — Confidence levels in (3 versus F for various dust approximations with the closed-box approximation. All the contours represent 
95% 2-para. confidence limits. The dotted line represents the result from marginalizing over attenuation = 0.2-0.55 and n = 0.8-1.2 
(Sec. 3.4); the upper solid line from the fiducial A-u = 0.35 and n = 1.0 (effective average over many galaxies); the dashed line from fixed 
Av = 0.75 and n = 0.65 (star-forming galaxies), and; the lower solid line from fixed Av = 0.5 with an MW extinction law (Pei 1992). 



Table 3 

Attenuation-corrected Ha luminosity densities 



reference 


log(LHa//iWMpc-^) 


notes 


Gallego et al. (1995) 
Trcssc & Maddox (1998) 
Glazebrook et al. (2003) 


32.39 

32.74 
32.77 


z « 0.03, UCM objective-prism survey, luminosity function 
z K, 0.20, CFRS, luminosity function 
2; « 0.03, SDSS, cosmic spectrum 



(Glazebrook ct al. 2003 included this uncertainty which 
amounted to 0.1-0.15 in the log result), (ii) uncertainties 
in the attenuation curve (e.g. Caplan & Deharveng 1986) 
which affect the conversion from a reddening measurement 
to absolute attenuation, and (iii) AGN contamination. 
Both the Gallego et al. (1995) and Tresse & Maddox (1998) 
results are based on obtaining the emission-line luminosity 
function from spectra. However, the selection criteria of 
the surveys are different, an emission-line objective-prism 
survey and /-band photometry, respectively. Glazebrook 
et al. (2003) took a different approach of summing SDSS 
spectra to form a cosmic optical spectrum before calcu- 
lating the Ha and H/3 luminosity densities. We take the 
mean of these three results for our fitting: 

log(LHa//iWMpc-^) =32.63 ±0.20. (6) 

This is appropriate since the Tresse^ & Maddox and Glaze- 
brook et al. results may be too high due to AGN con- 
tamination which are enhanced by aperture affects. The 
spectra taken through a fiber are normalized to a broad- 
band filter and since AGN are centrally concentrated, any 
emission line luminosity due to them will be over en- 
hanced. However, Glazebrook et al. did measure the weak 
Ol 6300 A line, suggesting that the AGN contribution was 
only a few percent at most. The Gallego et al. result is 
probably too low because of the small survey (i.e. large 
scale structure) and/or because it only includes EW>10 A 



(Ha, emission positive) galaxies. 

To compare this average measurement with model pre- 
dictions, we use the output from PEGASE for the Ha 

flux without dust attenuation."'' As with the earlier fit- 
ting, the synthetic spectra are normalized to the °-^r band 
so we are really comparing the ratio of Ha and ^'^r band 
fluxes between the models and the data. A correction is 
made for the fact that the measured "'^r band flux includes 
dust attenuation because we are fitting to an attenuation- 
corrected Ha measurement. We use the fiducial attenua- 
tion of Ay = 0.35. In other words, unlike the earlier fitting 
we now normalize the synthetic spectra to an attenuation- 
corrected luminosity density in the '^'^r band. 

Figure 13 shows the results from fitting to the average 
attenuation-corrected Ha luminosity density. The degen- 
eracy in P versus T is similar to that for the broadband 
luminosity densities. In addition, the best-fit IMF slope, 
r in the range 0.9 1.5 (68% conf. assuming a = and 
averaging over the metallicity approximations), is in good 
agreement with the those results. 

5. CONCLUSIONS 

In this paper, we present the results of fitting spectral 
synthesis models with varying IMFs to local luminosity 
densities. 

• A good fit is obtained to the measurements of Sul- 



® Some fraction of Lyman continuum photons are assumed to be absorbed by dust rather than gas according to the prescriptions of Spitzer 
(1978), normalized so that 30% of the photons are absorbed by dust at solaj: metallicity (Fioc & Rocca-Volmerange 1999). 
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r (IMF slope) 



Fig. 13. — Confidence levels in /3 versus F from fitting solely to the Ha luminosity density (Eqn. 6, relative to the " -'r band). The 
solid lines represent the contours for the closed-box approximation while the dotted lines represent the contours for the constant-metallicity 
approximation. The contours represent 68% 1-para. and 95% 2-paTa. confidence levels. 



livan et al. (2000), Blanton et al. (2003b) and Cole 

et al. (2001) if -band (compilation SBC^f) with a 
best-fit metallicity of around solar. If we fit to the 
measurements of the first two papers and Huang 
et al. (2003) (compilation SBH), the best fit metal- 
licity is greater than twice solar and the results are 
inconsistent with solar metallicity at the 99.7% con- 
fidence level. 

• The data can be well fit by a universal IMF and, 
therefore, there is no need to invoke IMF variations. 

However, this provides only a weak constraint on 
the invariance of the IMF because of significant de- 
generacies associated with this type of modeling. 

• The best-fit universal IMF slope marginalized over 

a significant range of cosmic SFH (0.5 < f3 < 4.0, 
-2 < a < 2) is r = 1.15 ± 0.2 based on an 
average between a closed-box approximation and 
a constant-metallicity approximation around solar 
(using compilation SBav). Our results are in good 
agreement with the Salpeter IMF slope. 

• A strong upper limit of F < 1.7 is obtained with 
99.7% or 95% confidence depending on the metallic- 
ity approximation. This rules out the Scalo (1986) 
IMF for a universal IMF since the mass fraction of 
stars above 10 Mq is similar to our parameteriza- 
tion with F = 1.8 (Fig. 1, see also Table 2). A sim- 
ilar conclusion was reached by Madau et al. (1998) 
from fitting to luminosity densities over a range of 
redshifts and by Kennicutt et al. (1994) from fit- 
ting to Ha EW-color relations for galaxies. Madau 
et al. found a Salpeter IMF or an IMF with F = 1.7 
provided adequate fits. Here, we find that the lat- 
ter slope does not provide a good fit to luminos- 
ity densities. This is principally because we have 
used more accurate local luminosity density mea- 
surements. If we increase our la uncertainties by 



0.05 at all wavelengths then we obtain F < 1.7 with 
80% confidence. 

• The stellar mass density of the universe is in the 
range 1.1-2.0 xlO"'^ based on marginalizing 
over cosmic SFH and our parameterization of the 
IMF. The current SFR density is in the range 0.7- 
4.1 xlO"2 /iMgyr^^Mpc"^. The total bolometric 
stellar emission (0.09-5 /um) is known more accu- 
rately, naturally because we observe light and not 
mass, and is in the range 1.2-1.7 xlO^^ /i WMpc"^ 
derived from the fitted PEGASE model spectra (the 
mass-to-light ratio is ~ 0.9 1.4 M0/L0). We find 
that our dust model can reproduce the estimated to- 
tal dust emission (0.3-0.7 xlO^^ /iWMpc"^, ~5- 
1000 /im) scaled from the far-IR luminosity density 
(Saunders et al. 1990) if we have A2000 < 1 (< 60%). 
Note that this represents a limit on the cosmic spec- 
trum attenuation, i.e., a luminosity- weighted aver- 
age by flux (not magnitudes). 

• Fitting to the local Ha luminosity density provides 
a similar result for the IMF slope (F = 1.2 ± 0.3). 
This provides some evidence that our upper mass 
cutoff of 120 M0 is a reasonable approximation be- 
cause the sensitivity of the Ha flux to massive stars 
is different to that of the mid-UV to optical fluxes. 
More accurate measurements could test this upper 
mass limit. 

• The quantitative results on F rely on the accuracy 
of the population synthesis model (PEGASE). An 
alternative, qualitative result would be that there 
is consistency between the theory of evolutionary 
population synthesis (evolutionary tracks and stel- 
lar spectra) and the measurements of luminosity 
densities, cosmic SFH and an average MW IMF de- 
rived from stellar counts (e.g. F « 1.3 from Kroupa 
2001). 
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Greater constraints can be placed on a universal IMF 
both by improved accuracy in local luminosity density 
measurements (in particular, z = 0.1 to match the multi- 
wavelength SDSS MGS) and by improved accuracy of di- 
rect measures of cosmic SFH with redshift (in particular, 
z = to z ~ 1-2). For the direct tracing of cosmic SFH, 
it is important that the UV is measured at the same rest- 
frame wavelength in order to avoid IMF dependency. In 
other words, we need an IMF-independent cosmic SFH in 



order for the local luminosity densities to accurately con- 
strain a universal IMF or the IMF-dependency should be 
quantified. 

We thank the referee Eric Bell for constructive sugges- 
tions, and also Timothy Heckman, Benne Holwerda, Rose- 
mary Wyse and Andrew Blain for helpful comments. We 
acknowledge generous funding from the David and Lucille 
Packard foundation. 
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